Comment on "Rovibrational quantum interferometers and gravitational waves 
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In a recent paper, Wicht, Lammerzahl, Lorek, and Dittus [Phys. Rev. A 78, 013610 (2008)] 
come to the conclusion that a molecular rotational- vibrational quantum interferometer may possess 
the sensitivity necessary to detect gravitational waves. We do not agree with their results and 
demonstrate here that the true sensitivity of such interferometer is many orders of magnitude worse 
than that claimed in the mentioned paper. In the present comment we estimate the expected energy 
shifts and derive equations of motion for a quantum symmetric top (diatomic molecule or deformed 
nucleus) in the field of gravitational wave, and then estimate the sensitivity of possible experiments. 



The direct detection of gravitational waves with an ex- 
perimental apparatus located on or near the Earth is one 
of the forefront goals of modern science. Several large- 
scale detector projects, which are presently being either 
built or designed, promise the first glimpses of this form 
of radiation that most certainly must permeate the fabric 
of the Universe. Given the scale of these experiments, it 
is reasonable to ask whether gravitational waves might 
be more simply detected. In a recent paper, Wicht et 
al. [H suggest that it is possible to detect such waves 
by measuring the energy shifts of rovibrational levels of a 
diatomic molecule. We show here that any such effect is 
many orders of magnitude smaller than predicted in [lj . 

Following Ref. [l[ we start with the problem of the 
hydrogen atom in the field of a weak gravitational wave 
propagating along the z axis, h^{t — z). The so-called 
TT gauge, where 



hu 



i-22, 



112 



i-21, 



hop 



0. 



(1) 



is convenient for the description of this wave. Wicht et 
al [l| base their analysis on the covariant Klein-Gordon 
equation. Expansion of this equation in powers of hij and 
1/c (c is speed of light) gives the following nonrelativistic 
hamiltonian for the electron of mass m and charge — e 



H 



2m 



(Sij - hij)didj - eA Q , 



(2) 



and where the field of a Coulomb center with charge e in 
the presence of the gravitational wave is 



A. = - 1 



^ijXi Xj 

~2r r ~ 



(3) 



This is the same hamiltonian used in [l[ , which we agree 
is valid. 
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Let us now calculate the energy shifts generated by the 
gravitational wave. The total system hamiltonian is 



H = 


H Q + 6H 






H = 


p 2 e 2 


2m r 






SH = 


2m J J 



~2~r~* 



(4) 



Naively, the energy shift is of the order of the perturba- 
tion 



Se ~ ( hijdidj) 



c f%ijXiXj 
2^3 



}-hRy, (5) 



where Ry stands for Rydberg, a typical atomic energy. 
However, we claim that the shift is much smaller, on the 
order of 



Se 



(hu) 2 



(6) 



where u> is the frequency of the gravitational wave, which 
is typically much lower than atomic or molecular frequen- 
cies of interest. 

To prove this we perform the following transformation 
of coordinates x — > y in the hamiltonian @ , (0]) • 



y, 



(7) 



This is the transformation into the frame freely falling 
in the gravitational wave. After this transformation, 
the interaction with the gravitational wave is completely 
removed from Eqs. (J2|) and (J3j> - However, this does 
not mean that the interaction with the gravitational 
wave vanishes completely. Consider the time-dependent 
Schrodinger equation 



ih 



d_ 

di 



yj(t,x) = Hip{t,x) . 



(8) 



The partial time derivative is taken at constant x. Under 
the transformation in Eq.((7|) the left hand side of Eq.® 
is changed as 
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d_ 

di 



dt 



dyj d 
dt dy l 



d_ 

dt 



2 



Therefore Schrodinger Eq. {HJ is transformed to 



H o(y) + -jhikVkPi 



ip(t,y), (10) 



where the equation is written in the y-coordinates and 
the time derivative is necessarily taken at y — constant. 
Hence, the new hamiltonian is 

H' = H + -hikVhPi = H + ^mh ik y k yi , (11) 
where we have used pi — miji . The new perturbation is 



SH' = ^mii ik y k yi = \^mh ik y k y, 
-> -^m'h lk y k y l . 



1 



-.•mKkVkVi 



(12) 



In this equation we have thrown away the full time 
derivative because it influences neither the equations of 
motion nor the physical results. Equation (|12p justi- 
fies our estimate (|6]) for the energy shift or for a typ- 
ical matrix element, if we substitute hi k ~ uj 2 h and 
yt,y k ~ og = h 2 /me 2 , where as is the Bohr radius. 

Obviously, the above statements are valid not only 
for a hydrogen atom. One can perform the transfor- 
mation 10 for any quantum system in a gravitational 
wave; furthermore one can certainly work in terms of 
original ^-coordinates as well. For example, for the prob- 
lem of the gravitational radiation of atomic hydrogen, 
the equivalence of the two approaches has been demon- 
strated explicitly in Ref. 0. However, in the case of 
low frequency gravitational waves, calculations in the 
^-coordinates ("x-frame" calculations) are "dangerous" 
because there is a dramatic exact cancelation (or com- 
pensation) between different large terms in the pertur- 
bation SH given by Eq. Q and approximations can 
easily destroy this exact compensation. Because the re- 
sult in [l[ does not contain the frequency of the gravita- 
tional wave and parametrically coincides with the naive 
estimate ([5]), we believe that this compensation was de- 
stroyed in Ref. [l[ due to approximations in the molecular 
calculations for the ionic molecule HD + . 

One can derive Eq. Q12p in a much simpler way using 
a frame freely falling in the gravitational wave from the 
very beginning. In this frame the center of mass of the 
system remains at rest. The motion of any particle of this 
system with respect to the center of mass is described by 
the well-known covariant equation (see, e.g., [1,01) 



Ds 2 



+ R% VT u p u T rf 



= 0. 



(13) 



In the general case, this equation for the so-called 
geodesic deviation rf" has the following meaning: Par- 
ticles a and b move along two close geodesies x£(s) and 



- and 



Xj(s); then, by definition, r? A '(s) 
?i M = dx 1 ^ (s) / ds . 

In our problem, let the index b label the coordinate of 
center of mass which is at rest in the chosen reference 



frame. Then — (1, 0, 0, 0), and in the weak field of the 

gravitational wave h mn = hmn exp[— iuj(t — z)] we arrive 
at the following equation of motion for particle a: 

1 



to h„ 



(14) 



In the hydrogen atom (or hydrogen-like ion), where 
the coordinates of electron, center of mass, and proton 
(or nucleus), x e , x c , and x p , respectively, lie on a straight 
line, we arrive immediately at the following equation of 
motion for the relative coordinate x = x e — x p : 



1 2, 

= LU h r , 



(15) 



The corresponding interaction hamiltonian for the hydro- 
gen atom is therefore 



Hint — ^ [1> ^ W h mn X m X n , 



(16) 



mP 



Let us note here that because h mn depend on t — z, but 
not on x m ^ n , the relations (|15p and (|16[) are in exact 
correspondence. Eq. (p~6|) coincides with Eq. (fT2| . 

The molecular ion HD + considered in Ref. [l(] con- 
sists of three constituents: an electron, a proton, and a 
deuteron. To simplify the description, we note that, with 
m e m p , m d , one can assume that the coordinates of 
proton, center of mass, and deuteron, x p , x c , and x d , re- 
spectively, lie on a straight line, so that for the relative 
coordinate x = x p — x d we have the same equation, (|15l) . 



Therefore, the interaction hamiltonian for the HD + ion 
can be written as 



Hi, 



^ m i l -E m -E n ) 



(17) 



-.d p 



mP 



Tfl . 



We have neglected here the formally-arising term 



1 e 2 

— m ur 
4 



because the assumption m e -c mP has been made already 
for the derivation of hamiltonian (|17[) : furthermore, the 
contribution of this term to the discussed effect is rela- 
tively suppressed as m e /m p . 

The typical frequency shift for a HD + ion, as estimated 
with hamiltonian in Eq. (|17p . is 



1 



5v ~ —m p LU 2 ha 2 n 



(18) 



This estimate is valid as long as the frequency to of the 
gravitational wave is much less than that of the molecular 



3 



vibrations and rotations. For the dimensionless ampli- 
tude h ~ 10 -19 and frequency uj ~ 10 4 Hz, we estimate a 
frequency shift as ~ 10~ 23 Hz. This estimate is well be- 
low the corresponding estimate ~3x 1CU 5 Hz presented 
in Ref. [l[ . Moreover, this frequency shift varies with time 
together with hit — z). With this time-dependence of the 
frequency shift, the observation of the effect seems even 
less realistic because high precision spin precession mea- 
surements generally require considerable averaging time 
(days to years). 

Let us make some additional remarks about the be- 
havior of a molecule (or a deformed nucleus) under the 
influence of the hamiltonian (|17p . This hamiltonian in- 
fluences both the vibrational and rotational degrees of 
freedom of the diatomic molecule. The classical (non- 
quantum) vibrational dynamics in the gravitational wave 
are discussed in several textbooks, see e.g. Ref. and 
quantum mechanical treatment does not bring significant 
differences in the dynamics. Here we derive the implica- 
tions of Eq. (jTTJ) for rotational dynamics. The classical 
equation of motion for angular momentum L = r x p 
follows from Eq. (JTTJ) 



Li 



1 



-L0 2 €iklhi n lk, 



(19) 



where is is the tensor of inertia of the molecule. Only 
the traceless part of Ik n contributes to Eq. I|19p. and 
the molecule is a symmetric top. As usual we denote 
by L | the moment of inertia for rotations around axis of 
the top, and by I± the moment of rotations around a 
perpendicular axis. 1 Then 



h 



1 



(2/i 



n k n Tl 



(20) 

Clearly Eqs. (fT9|) and (|20|) are valid for both the classical 
and quantum cases. Let us consider a molecule in a ro- 
tational quantum state with a given angular momentum 
J, L — hJ, J 2 = J(J + 1). Then (see e.g. @) 



' An 



I I JkJn + Jn-h ~ ^ J (J 



I = Iu-Ij 



l)4ri 

o 

3n 2 - j(j + i) 



- . (21) 

J{2J- 1)(J+ l)(2J + 3) V ; 

Here we keep only the traceless part of the inertia tensor, 
and f2 is projection of J on axis of the top. Hence Eq. 
(Ti"9l) is transformed to 



(22) 



X 2 

•h = tiklhln {JkJn + JnJk} 

The corresponding effective Hamiltonian is 

1 
2 



Heff — —-^uhlnJlJn 



(23) 



1 For a molecule or deformed nucleus one must set /| | =0. Never- 
theless, we keep /j i to stress that only traceless part of the inertia 
tensor contributes to the spin precession. 



Eqs. (|21j) and ([23)1 are equally applicable to a diatomic 
molecule or to a deformed nucleus, and describe transi- 
tions or energy shifts between magnetic sublevels. For 
the ground state J = fl = S, S is spin of the system. We 
assume that S > 1/2. Hence Eq. (|22"|) is transformed to 
equation of spin precession in the field of gravitational 
wave 



Si 



2h(S + 1)(2S + 3) 



(S k S n + S n S k ) . (24) 



This is a general equation for a free spin. 

Given that a magnetic resonance frequency shift due 
to (|2"4"|) is extremely small, we alternatively consider the 
case of a solid, where there is a possibility of a much 
larger nuclear electric quadrupole field shift. This shift 
is induced by the deformation of the lattice due to the 
gravitational wave. For a sample with size ~ 1 m and 
for uj > 10 4 Hz the frequency of the gravitational wave 
disappears from the equations because the frequency of 
the lowest acoustic standing wave ujq is lower than to. 
The cancellation of the w-dependence is similar to that 
for the LIGO experiment discussed below, see eqs. (|2T| . 
(|2"5)l . Hence the maximum possible nuclear quadrupole 
resonance frequency shift is 



S - 10 9 /i Hz 



(25) 



for the case of deformed nuclei. 

Unfortunately the frequency shift 10 -10 Hz given by 
Eq. (|2"5j) with h ~ 1CP 19 is still very small. The current 
state of the art nuclear spin precession frequency mea- 
surements accuracy (for a gas of 199 Hg atoms used in a 
permanent electric dipole moment (EDM) search [7|) is 
barely at this level after several years of measurement 
and averaging, and this experimental accuracy far ex- 
ceeds that of any other experiment performed to date. 
Of course, the possibility of a macroscopic sample with 
a large number of participating spins suggests that some 
increase in sensitivity might be possible, but there is a 
considerable loss in sensitivity due to large-sample mea- 
surement limitations. The only plausible measurement 
technique is the much less efficient monitoring of the sam- 
ple's precessing magnetization with a SQUID or other 
magnetometer, compared to the direct optical detection 
of the spin state that can be used with gaseous samples. 

Let us add a brief remark about gravity wave detectors 
based on strain measurements (such as LIGO). These are 
also described by Eq. ([14]) . Suppose a mass (such as 
an interferometer mirror) is positioned in equilibrium at 
x = y = z = 0, and another mass (the second mirror) 
is positioned in equilibrium ai x — L, y = z = 0. The 
mass at the origin is unaffected by the gravity wave, the 
response of the mass at x = L is given by 



x + oJq(x — L) 



1 



-oj 2 h+(t)L, 



(26) 



where ujq is the natural oscillation frequency around the 
equilibrium point, h + = hn = — /122, and we neglect 



4 



damping. For a wave h + (t) — h+e lult the solution is 



x = L + - 



1 CO 2 



h^Le 



2uj 2 -uj 2 + 



(27) 



LIGO is a "free-mass" detector, which means that ujq <c 
lu (uj ps 4 Hz for Advanced LIGO [1]). This gives the 
strain 



AL/L = -h+ e 



2 



(28) 



The strain in oscillator modes with frequencies much 
greater than the gravitational wave frequency is sup- 
pressed by uj 2 /iv 2 . Indeed, the more rigid the system, 



the less sensitive it is to an external force. In a molecule, 
for example, lu is on the order of THz, which gives an 
enormous strain suppression for gravitational wave fre- 
quencies io <J 10 4 Hz. 



Acknowledgements 

The authors would like to acknowledge useful discus- 
sions with David DeMille, Dmitry Budker, and Max 
Zolotorev. The work by LB. Khriplovich was supported 
by the Russian Foundation for Basic Research through 
Grant No. 08-02-00960-a. 



[1] A. Wicht, C. Lammerzahl, D. Lorek, and H. Dittus, Phys. 
Rev. A 78, 013610 (2008). 

[2] S. Boughn and T. Rothman, Class. Quantum Grav. 23, 
5839 (2006); lgr-qc/0605052"1 v.2 (see, in particular, sec- 
tion 5 therein)! 

[3] J. Weber, General Relativity and Gravitational Waves, In- 
terscience Publishers Inc., New Jork, 1961. 

[4] A. P. Lightman, W.H. Press, R.H. Price, S.A. Teukol- 
sky, Problem Book in Relativity and Gravitation, Prince- 
ton University Press, Princeton, New Jersey, 1975, prob- 
lem 9.23. 



[5] C. W. Misner, K. S. Thorne, J. A. Wheeler, Gravitation, 

San Francisco : W. H. Freeman, 1973, section 37.5. 
[6] L.D. Landau and E.M. Lifshitz, Quantum mechanics: non- 

relativistic theory, Oxford, New York, Pergamon Press, 

1965, section 119, problem 1. 
[7] W. C. Griffith, M. D. Swallows, T. H. Loftus, M. V. Ro- 

malis, B. R. Heckel, and E. N. Fortson, Phys. Rev. Lett. 

102, 101601 (2009). 
[8] Advanced LIGO Reference Design Document, LIGO 

M060056-08-M, 2007. 



